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Abstract 

The cascading gauge theory of Klebanov et. al reahzes a soluble example of gauge/ string 
correspondence in a non-conformal setting. Such a gauge theory has a strong coupling 
scale A, below which it confines with a chiral symmetry breaking. A holographic 
description of a strongly coupled cascading gauge theory plasma is represented by a 
black brane solution of type IIB supergravity on a conifold with fluxes. A characteristic 
parameter controlling the high temperature expansion of such plasma is Q ~ (in ^) ^ . 
In this paper we study the speed of sound and the bulk viscosity of the cascading 
gauge theory plasma to order 0(0^). We find that the bulk viscosity satisfies the 
bound conjectured in arXiv:0708.3459. We comment on difficulties of computing the 
transport coefficients to all orders in Q. Previously, it was shown that a cascading 
gauge theory plasma undergoes a first-order deconfinement transition with unbroken 
chiral symmetry at Tcruicai = 0.6141111(3)A. We show here that a deconfined chirally 
symmetric phase becomes perturbatively unstable at TunstaUe = 0.87 49 {0)TcriUcai- Near 
the unstable point the specific heat diverges as cy ~ |1 — 
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1 Introduction 

Consider J\f — 1 four-dimensional supersymmetric SU{K + P) x SU{K) gauge theory 
with two chiral superfields Ai,A2 in the {K + P, K) representation, and two fields 
Bi, B2 in the {K + P,K). This gauge theory has two gauge couplings gi, g2 associated 
with two gauge group factors, and a quartic superpotential 

W ^Tj:{ABjAkBe)e%^'. (1.1) 

When P = above theory flows in the infrared to a superconformal fixed point, 
commonly referred to as Klebanov-Witten (KW) theory [1]. At the IR fixed point 
KW gauge theory is strongly coupled — the superconformal symmetry together with 
SU (2) X SU (2) X f/(l) global symmetry of the theory implies that anomalous dimensions 
of chiral superfields 7(Aj) = 7(-Bi) = — |, i.e., non-pcrturbatively large. 

When P 7^ 0, conformal invariance of the above SU {K + P) x SU (K) gauge theory 
is broken. It is useful to consider an effective description of this theory at energy scale 
H with perturbative couplings gi{fi) <^ 1. It is straightforward to evaluate NSVZ beta- 
functions for the gauge couplings. One finds that while the sum of the gauge couplings 
does not run 

d (. 0, (1.2) 



the difference between the two couplings is 

where A is the strong coupling scale of the theory and is an anomalous dimension 
of operators TrAiBj. Given (1.3) and (1.2) it is clear that the effective weakly coupled 
description of SU {K + P) x SU {K) gauge theory can be valid only in a finite- width 
energy band centered about ^ scale. Indeed, extending effective description both to 
the UV and to the IR one necessarily encounters strong coupling in one or the other 
gauge group factor. As beautifully explained in [2], to extend the theory past the 
strongly coupled region(s) one must perform a Seiberg duality [3]. Turns out, in this 
gauge theory, a Seiberg duality transformation is a self-similarity transformation of the 
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effective description so that K ^ K — P a,s one flows to the IR, oi K ^ K + P a,s the 
energy increases. Thus, extension of the effective SU {K + P) x SU (K) description to 
all energy scales involves and infinite sequence - a cascade - of Seiberg dualities where 
the rank of the gauge group is not constant along RG flow, but changes with energy 
according to [4-6] 

K = K{f^)^2Pnni^, (1.4) 

at least as » A. To see (1.4), note that the rank changes by AK ~ P as PA (in ^) ~ 
1. Although there are infinitely many duahty cascade steps in the UV, there is only a 
finite number of duality transformations as one flows to the IR (from a given scale /i) . 
The space of vacua of a generic cascading gauge theory was studied in details in [7] . In 
the simplest case, when K{ii) is an integer multiple of P, the cascading gauge theory 
confines in the infrared with a spontaneous breaking of the chiral symmetry [2] . 

Effective description of the cascading gauge theory in the UV suggests that it must 
be ultimately defined as a theory with an infinite number of degree of freedom. If so, 
an immediate concern is whether such a theory is renormalizable as a four dimensional 
quantum field theory, i.e., whether a definite prescription can be made for the com- 
putation of all gauge invariant correlation functions in the theory. As was pointed out 
in [2], whenever gsK{ij) ^ 1, the cascading gauge theory allows for a dual holographic 
description [8, 9] as type IIB supergravity on a warped deformed conifold with fluxes. 
The duality is always vahd in the UV of the cascading gauge theory; if, in addition, 
QsP ^ 1 the holographic correspondence is valid in the IR as well. It was shown 
in [10] that a cascading gauge theory defined by its holographic dual as an RG flow 
of type IIB supergravity on a warped deformed conifold with fluxes is holographically 
renormalizable as a four dimensional quantum field theory. 

Cascading gauge theories provide a soluble realization of the holographic gauge 
theory/string theory duality in non-conformal setting. A way to construct four dimen- 
sional examples of non-conformal gauge theory /string theory correspondence is to start 
with an AdS^/CFT^ duahty and to deform it by relevant operators of the CFT4. An ex- 
ample of such construction is the gauge/string duahty for the J\f — 2* supersymmetric 
gauge theory [11-13]. On the contrary, the scale in a cascading gauge theory is intro- 
duced via a dimensional transmutation of the gauge couplings (1.3). We would like to 
understand in details the hydrodynamic properties of strongly coupled non-conformal 
gauge theory plasmas [14]. A lot is known about thermodynamics/hydrodynamics 
of strongly coupled mass deformed conformal gauge theories from the perspective of 
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gauge/string correspondence [15-19]. The thermodynamic/hydrodynamic analysis of 
the cascading gauge theory plasma are substantially more difficult. The equilibrium 
thermodynamics of the cascading gauge theory plasma in the deconfincd chirally sym- 
metric phase is well understood by now [20-22]. Since at zero temperature a cascading 
gauge theory confines with a chiral symmetry breaking, it is conceivable that there 
is a finite-temperature deconfined phase of the theory, with broken chiral symmetry^ . 
Whether or not such a phase exists is an open question [24]. In case of hydrody- 
namic transport coefficients^, the shear viscosity was shown to satisfy the universal 
bound [26-28], and the bulk viscosity was computed to leading order at high temper- 
ature [29]. 

In this paper we study propagation of sound waves in the strongly coupled de- 
confined chirally symmetric phase of the cascading gauge theory plasma. In the dual 
gravitational description this involves computation of the dispersion relation of the 
lowest quasinormal mode in the sound channel [30] of the black hole solution numeri- 
cally constructed in [22]. Ideally, we would like to do the analysis at any temperatures 
(at least above the dcconfinement transition), much like it was done for the M = 2* 
plasma in [19]. Unfortunately, in section 3 we show that technical difficulties in the 
present framework does not allow us to achieve that goal. Thus, we resort to pertur- 
bative high temperature computations. The small parameter of the high temperature 
expansion is ^ [21], where ~ is roughly^ the effective rank of the cascading 

plasma (1.4) at the IR cutoff scale, set by the temperature. We compute both the 
speed of the sound waves and the bulk viscosity of the cascading gauge theory plasma 
to order O {j^- Previously such analysis were done to order O (y^^ [29]. At the 
order reported, the bulk viscosity was shown to saturate the bound proposed in [31]. 
Higher order corrections to the bulk viscosity of the cascading plasma presented here 
show that the bound [31] is satisfied. Alternative way to compute the speed of sound 
is to use the equilibrium equation of state 

Using (1.5) and shghtly extending the computations in [22], we can evaluate the speed 
of sound to temperatures down to the dcconfinement transition and below. The com- 



^Such a phase was observed in 4 -|- 1 dimensional supersymmetric SU{Nc) gauge theory with 
fundamental quarks compactified on a circle [23]. 

^Hydrodynamics of closely related models was recently discussed in [25]. 
"^A precise definition of iC* is given below. 
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parison between the perturbative high temperature analysis and the exact one indi- 
cates that the former is convergent for ^ ~ 2 ■ ■ - 3, correspondingly to temperatures 
T ~ (1 • • • 1.5) A — which is about twice the critical temperature of the deconfinement 
phase transition Tcriucai = 0.6141111(3)A [22]. Thus, although we find a relatively small 
bulk viscosity in the high temperature regime ^ ^ we can not (reliably) evaluate the 
bulk viscosity in the vicinity of the deconfinement transition. It it clear though that 
since the deconfinement phase transition is of the first-order (in the 't Hooft limit), the 
bulk viscosity will remain finite at the transition point [32]. Finally, we find that chirally 
symmetric phase becomes perturbatively unstable at T = Tunstabie = 0.87487(7)Tcriticai 
— exactly at this temperature vanishes, and extending this phase to lower values of 

leads to < 0. The critical behavior at the unstable point in the cascading plasma 
is identical to the one found m. M — 2* plasma with mass deformation parameters 
nif < rrib [19,31]. We comment more on the instability in section 5. 

The technical aspects of the computations are presented in section 2. The reader 
interested only in the results should consult section 4. 



2 Supergravity dual to deconfined cascading plasma 

The holographic dual to a deconfined chirally symmetric phase of a cascading gauge 
theory plasma at equilibrium is given by a black hole solution in a singular Klebanov- 
Tseytlin (KT) geometry [33]. It has been discussed previously in [4,20-22]. We follow 
the notations of [22]. A black hole metric^ : 



dsl^ =h-^''^{2x - x^Y^I'^ (-(1 - xfdt^ + dx\ + dxl + dx^^ + G^^{dxf 

a=l 

where /i, ji and /s are some functions of the radial coordinate x. There is also a dilaton 
g{x)^ and form fields given by 

F3 = P A (e^i A e^^ - e^^ A e^J , P2 = ^ (e^i A e^^ - ee^ A e^^) , 

(2.2) 

-^5 = -^5 + *-^5 , -^5 = A egj A e^^ A A e^^ , 



*The frames {eg^^e^^ are defined as in [10], sucli tliat tfie metric on a unit size T^'^ is given by 



(4)+eLi( 
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where K is a function of the radial coordinate x. Without loss of generality we can 
set asymptotic string coupling to one. We use the following parametrization for the 
solution in perturbation theory in ^ : 

- § + Ig { (e-(-) - lv.nix)) } , (2.3) 



f2ix) =do + ao ^1 (^^^ (^-2^2nix) + ri2nix) + ^\2nix)^ | , (2.4) 

fsix) = So + ao { (^-'2^2n{x) + ri2n{x) - ^A2„(x)^ I , (2.5) 

K{x) -K. + K^f^!^ K2n{x) I , (2.6) 

9i^) = ^ + f2[[^^ C2nix)Y (2.7) 

The advantage of this parametrization is that the equations for {^2n, V2n, ^2n, C2n} de- 
couple, once the (decoupled) equation for K2n is solved, at each order n in perturbation 
theory. 

Gauge invariant fluctuations 

{Zh , Zf , , , Zk} 
of the background metric (2.1) and the scalar fields 

|/ = iln/i + ^ln/3 + ^ln/2, u; = lln|, $ = ln^, k| 

of the effective five-dimensional gravitational description of the sound channel quasi- 
normal modes were studied in details in [29]^. The incoming wave boundary conditions 
on all physical modes imply that 

Zh{x) = (1 - xY'^'zHix) , Zf{x) = (1 - x)-'^Zf{x) , 

Z^{x) = (1 - xy'^z^ix) , Z^{x) = (1 - xy'^z^ix) , (2.8) 

Zk{x) = (1 - xY'^^ZKix) , 



^See eq. (63) of [29] for the definition of the gauge invariant fluctuations and eqs. (64)-(68) of [29] 
for the corresponding equations of motion. 
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where {zh, Zf, z^, z^, zk} are regular at the horizon; we further introduced 



W = 



00 



q = 



27rr' 27rr' ^^'^^ 

T is the equihbrium temperature of the plasma, and {u!,q — \q\} are the frequency and 
the momentum of the sound quasinormal mode. There is a single integration constant 
for these physical modes, namely, the overall scale. Without the loss of generahty the 
latter can be fixed as 



zh{x) 



(2.10) 



In this case, the sound dispersion relation is simply determined as 



zh{x) 



= 0. 



(2.11) 



X~^0-L 



The other boundary conditions (besides regularity at the horizon and (2.11)) are 





= , z^{x) 


— , z^{x) 


= 0, zxix) 











= 0. 



(2.12) 



Let's introduce 



Zh = 



2\ " 



?2 \ " 



p2\ " 



.^g{(S"4:2}^M|{(a"4:4, 

n=0 *^ > n=0 

°°j'/p2\" 1 °° 

*=^.f:l(0"4"i}+M^.E{( 



J") 



*,1 r ' 



(2.13) 



p2 



n=0 \ ■> n=0 

where the lower index refers to either the leading, oc q°, or to the next-to-leading, 
oc q^, order in the hydrodynamic approximation, and the upper index keeps track of 
the ^ high temperature expansion parameter. Additionally, we find it convenient to 
parametrize 



ASKS'-- 



n=0 
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where /3i „, /?2,„ are constants which are to be determined from the pole dispersion 
relation (2.11) 



(n) 
%,0 



= 0. (2.15) 



2.1 Equilibrium thermodynamics to order O (^^^ 

Cascading gauge theory thermodynamics for generic ^ was studied in [22]. Since we 
are able to compute the O (q^) sound wave dispersion relation only pcrturbatively in 

we need to construct the black hole backgrounds in the high temperature expansion 

as well. To order O (^^j this was done in [21], and extended to order O in [22]. 
Here, we extend the analysis to n = 4 in (2.3)-(2.7). 

Equations of motion for {K2n , ^2n , V2n , ^2n , C2n} take form 

o = ^'L + :^ + Jl':\ (2.16) 

Q-r," + _ ^V2n _ 2 , , _ 8^2n ^[2n] 

:r2(x-2)2 5"^''^" ?,x^{x-2Y^^'^^ ' ^^-^^^ 

(3x^-6x + 4)e^„ _ 2 pn] . . 

"-^2n+ ^(^_1)(^_2) 3'^2'^2n + A« , (2.18) 

n-A" + ^^2n g / / ^[2n] /2iq>, 

U-A2„+_^_^ x2(x-2)2 ^^^2 /^Sn + ->'&,A ' 1^--^^^ 



= C + :;%r + 24 4„ + J-J?' , (2.20) 

where the source terms {J^^\ J^^\ J^^\ jfx\ J^^^} are functionals of the lower or- 
der solutions: K2m, C2to, ''?2to, A2m, C2m, with m < u. Explicit expressions for the source 
term functionals are available from the author upon request. The perturbative solu- 
tions to (2.16)-(2.20) must be regular at the horizon, and must have the appropriate 
KT asymptotics near the boundary. Beyond n — 1, these equations must be solved 
numerically. We apply the numerical strategy developed in [22]: 

■ Generically, the differential equations will have non-normalizable modes near the 
boundary^ x ^ 0+, and can generate singular Schwarzschild horizon as a; ^ 1_. Thus, 



^These modes are singular in case they are dual to operators of dimension larger than four — as 
for {f,Lo}; in other cases they modify the KT asymptotics, i.e., the parameters (such as a strong 
coupling scale) of the dual plasma. 
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we specify boundary conditions as a series expansion near the boundary and the hori- 
zon which exphcitly contain only normahzable modes. 

■ The total number of integration constants near the boundary and the horizon ap- 
pearing in normahzable modes precisely equals the total order of the system of ODE's. 
As a result the boundary value problem is well posed. 

■ We solve the resulting boundary value problem as detailed in section 5.2 of [22]. 

In what follows we present the horizon and the boundary expansion of the normal- 
izable modes for {K2n , ^2n , V2n , ^2n , (211} with n = {1, 2, 3, 4}. In numerical analysis 
we used expansion to order 0{x^^'^) (up to powers of \nx) near the boundary; and to 
order 0{y^^), y = 1 — x, near the horizon. Below, however, we present expansion at 
most to orders 0{x'^) and 0{y'^) — just what is enough to exhibit the dependence on 
all the integration constants. We solve the boundary value problem on the interval 
X e [5x-i 1 — 5x] with 5x = 10~^; we verified that our final results are insensitive to the 
precise choice of 5^, provided it is sufficiently small. 

Notice from (2.18) that ^2n always has a zero mode. Such a zero mode simply 
rescales (perturbatively in ^) oq, and has no effect on physical quantities [22]. In 
what follows we conveniently set this mode to zero near the boundary'' — for further 
details see (2.31), (2.42), (2.53) below. Naively, from (2.20), ^2,1 also always has a zero 
mode. The latter however is fixed by our choice of the asymptotic string coupling (2.7). 
Similarly, the zero mode of K^m see (2.16), modifies the strong coupling scale A of the 
cascading plasma. 

2.1.1 Order n = 1 
We find: 

~].n{2x-x^), (2.21) 

e2 = ^ln(2x-x2). (2.22) 

Even though it is possible to write down explicit analytic expressions for {772, A2, C2}, 
such expressions involve complicated polylogarithm functions, which slows down sub- 
sequent numerical computations. Thus, we opt to treat these fields numerically. 

^Of course, it is inconsistent to require the vanishing of the zero mode both near the boundary and 
the horizon. 
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Near the boundary, x — > 0+, we have: 

ri2 = - ^ + ^ln2 + ^Inx - + (^rjt'^ + ^Inx^ + 0{x^ Inx) , (2.23) 

X, ^^x + Xl'\'/' + 0{x') , (2.24) 

(2 =x (^Ci'° + ^lnx^ + 0{x^ Inx) . (2.25) 
Near the horizon, y — > 0+, we have: 

m ^vIh + 0(y') , (2.26) 

A2 + 0{y') , (2.27) 

C2 =Ci% + Oiy') . (2.28) 
Altogether at this order we have 6 integration constants 

K'°,Af ,C,^'°,<,,A°,,Cy, (2.29) 
which is precisely what is needed to specify a unique solution for {772, A2, C2}- 

2.1.2 Order 2 

Near the boundary, x ^ 0+, we have: 

(^4'^ - ^ In + 0(^^/2) , (2.30) 

U=^\nx + 0{x \nx), (2.31) 
00 

riA — 1 ln2H Inx + x ( Ci^'" H In 2 ° H In a; 

12 18 18 V 30^^ 360 60 30 ^ 12 

A?V/^ + (xJ" + (-5"" + ll^f ° - - S ^ - 1^^^^°) '"^ 

(2.32) 

A4 =x (^Kf - y + ^Cf + ^ In 2) + x^/^ Af + 0{x' In x) , (2.33) 

11 
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C4 -X (^Cr + (/^f - ^ + il^2 + ic^°) Inx) + 0{x' ]n'x) . (2.34) 
Near the horizon, y ^ 0+, we have: 



«4 ^Kl„ + 0{y^), 


(2.35) 


U HIh + + ^(^') ' 


(2.36) 




(2.37) 


A4 , + 0{y') , 


(2.38) 


C4 =C2% + Oiy') . 


(2.39) 



Altogether at this order we have 10 integration constants 

{l^2 jV2 X2 ) 5 'C2,/i 5 C2,/n ''?2,/i 5 ^2,h ! C2,/t} ) (2.40) 

which is precisely what is needed to specify a unique solution for {ac4, ^4, 774, A4, ^4}. 

2.1.3 Order n^?, 

Near the boundary, x — > 0+, we have: 

«6 -X + (-4° + ^ - i - ) Inx) + 0{x''') , (2.41) 

1,2 17 , 49 / 1 , 1 \ , / 47 1 no, 

776 = ln2H In 2 h ln2H Inx + x Ci In 2 

54 216 648 V 54 24^ V 1080 60^^ 

_ 1^2,0 In 2 + In 2 + — - -K^ - -C2' - - — In^ 2 

30 ^ 2160 360^ 60 ^ 30^ 30 ^ 120 

+ f 1^2,0 + 1 - 1^ Inx^ + (-—>cf - —X'A 

\% ^ 12^^ 24 9/ y V 225 ^ 225 ^ / 

I 2/ 4,0 , / 1^2,0 2,0 I 1 / 2,0\2 1 2,0i o I 89 4,0 19 i o ^2,0 

+ x 77o' + -Ci H U«9 ) K"> ln2H 77i' ln2 Ci 

V V9 ^ 15^ ^ ^ 45 ^ 72'^ 1440 1620^^ 

_ 11/^2,0 25121 _ 1^ 4,0^^2 - Irj^'O + —{C^'^Y + —C^'^ + —k'^ 
90 ^ 259200 18 3 540^^^ ^ 18^^ 15 ^ 

- — In^ 2 Vn x + i--Kf + -77^ - -C^° - — In 2 + ^\ In^ x 
180 J V 15 12 30^ 720 8640 / 

+ Y^ln'^)+C(^'/' Inx), 

(2.43) 
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Ae =x (-| In 2 + ^ + ^C^^ In 2 + In 2 + i In^ 2 - ^-^Cl' - f «f + ^Cf 
+ ^4°) +x^/^Xf + 0{x^ Inx), 



(2.44) 



-2,0 , / 5 2,0 5 2,0 , 1 /•2,0 , 2,0, 25 5 r, l,2r),l ^2,0 ■ 



Ce Cs' + --Cr - -1^2 + -C2' + i^T + — - - In 2 + - In^ 2 + -(r" In 2 

Sb I S3 gSl 3 2 3 -^g g g ^Sl 

-'- 2,0 1 o 1 1„ ™ 1 , /'/^/™2 1 2 



+ -Krin2j lna;J + In^a;) . 

(2.45) 

Near the horizon, y — > 0+, we have: 



«6 =4, + C»(y'), 


(2.46) 




(2.47) 


^6 + 0{y^) , 


(2.48) 


Ae =A° , + 0{y') , 


(2.49) 


Ce =C3% + 0{y') . 


(2.50) 



Altogether at this order we have 10 arbitrary integration constants 

which is precisely what is needed to specify a unique solution for {kq, ^q, rjQ, Ag, Ce}- 
2.1.4 Order 4 

Ultimately, to evaluate the speed of sound and the bulk viscosity at this order we will 
need only Kg and ,^8 solutions. 

Near the boundary, x 0+, we have: 

«8 -X + (^C^° + y - - - I + ^ In 2 - 1 In^ 2 - lc^° In 2 
-^4° In 2^ Inx^ +0{x'/^), 

(2.52) 
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n 


1 
1 


9 
Z 


Q 
O 


A 


9 




0.73675974 


-0.62226255 


-0.03784377 


A n 
In 


-0.01717287 


0.00534036 


-0.01064222 




A ' 


-0.87235794 


-1.11562943 


-1 c\ f\ r\ r\ c~\ ft c\ ft 

1.39008636 




Sn 


-0.15342641 


0.62226267 


-0.32514260 









0.62226259 


f\ A f\ f \ fi -t A ft ■% 

-0.42061461 


f\ f\ f\ f\ -t ft ft ct -t 

0.00816831 


j-Q 




-0.07981931 


f\ f \ 't ft ft -1 f\ 

0.01661150 


f\ f\ f\ fy f\ f\ Ck ^ fy 

-0.00920379 


^0 
^n,h 




0.01919989 


-0.05277626 


0.01385333 


'ln,h 


-0.14891337 


-0.21809464 


0.00213345 






0.16806881 


-0.14619173 


0.01639579 






-0.41123352 


0.33024116 


-0.07445122 





Table 1: Coefficients of the normalizable modes of the background geometry. See 
(2.29), (2.40), (2.51) and (2.56). 

/ 1 23 41 \ 

^» = (324 2 - 12% 2 + 1944 j ^ + ^) ■ 

Near the horizon, y ^ 0+, we have: 

Ks ^kI„ + 0{y') , (2.54) 

=e,,H + a,Hy' + 0{y') . (2.55) 
Altogether at this order we have 4 integration constants 

{K'A^ , K'4^h ■> ^4,,h ■> ^l,h} ■> (2.56) 

which is precisely what is needed to specify a unique solution for {ft^s^Cs}- 
2.1.5 Integration constants for the normalizable modes 

Here we tabulate (see table 1) the integration constants for the normalizable modes 
of t^2n, (,2n, V2n, ^2n, (2n with u = {1,2,3,4} obtained from solving the corresponding 
boundary value problems. 
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2.1.6 V , E and cj from equilibrium thermodynamics 
Using the results of [22] we can compute 



= ^(t^-«2,o1, |. = 7(1 + ^«2,o) , (2.57) 



sT 7 \12 "'7 ' sT 4\ 7 

where s is the entropy density^. Furthermore, the perturbative high temperature ex- 
pansion for 02,0 is given by^ 

7 P2 /7 35 7 7 , \ P4 / 35, ^ 35 20 7,2^ 

+ IZ^ + 1(2,0 ^ 2 + l«,2,o ^ 2 + 1(2,0 7 2,0 _ 35 2,o\ + ("ITS 
108 24^^ 12 ^ 12^^ 6 ^ 36^^ ) Kl \72 

I 7 2,0 , 175 2,0 35 2,0 1 o 35 2,0 1 r, ^75 175 2,0 35 20 

H G H in 2 Ci m2 Ci in 2 

12^^ 54 ^ 18 ^ 36^^ 324 108^^ 48 

+ 1 ln3 2 + 1(2'0 ln2 2 + -Kf In 2 + -k^ In^ 2 + 1(2^'° In 2 + -Acf 
96 48^^ 12 ^ 24 ^ 24^^ 6 ^ 



36^^ 18 y V i^l 



(2.58) 

Note that the coefficient of ^ must vanish [22] — numerically we find that it is 
oc 2 X 10-1°. 

The precise temperature dependence of K^, was determined in 



1, /647r4 sT\ 
^ = -ln X — . 2.59 

P2 2 V 81 ^ / 

Using (2.59) and the expressions for the pressure and the energy density from (2.57) 
we find 

2_g^_ 1 7 - 12a2,o - 6P2 ^ 

88 3 7 + 4a2,o + 2P2^ ' ^'"''^ 

Thus, given the perturbative high temperature expansion for 02,0 we can evaluate from 



^Note that expressions in (2.57) are valid for any temperature. 

^In order to evaluate the coefficient of ^ term we need the boundary expansions for 7^4 and A4 
to order 0{x). These expansions do not depend on the coefBcients of the corresponding normalizable 
modes {??4'° , A4'°}. 
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(2.60) the perturbative high temperature expansion for c 

4p2 



3cf = = 1 



+ 



3K 
10, . 1 

3 3 



3 3 3 ^3^^ 



^\(-8-^C^°ln2 



I 1 o ^ 1 2 o ^ 2,0 , 40 2,0 

H m2 m 2 k,' H Ci 

Q o o o g 



^2,0, o I 2,0 

-K2 m2 + -^'^2 



3^2 



+ (-i In^ - \>^r In 2 - In^ 2 - In 2 + ^(4°)^ + y /^^'°C^° 
+ ^(C^°)^ - ^ In^ 2 - 12 In 2 + ^C^° In 2 + 1^ + ^^^f In 2 + ^ In^ 2 



212 2,0 106^2,0^46^2,0 



92 



- — «o - — Ci' + yC2' + y«^3 



2,0 



3C3 



2,0 



2,0 



p8 



)10 



X5 



(2.61) 



Notice that (2.61) provides predictions for of (2.14) with n = 1, • • - 4. 

2.2 Speed of sound waves in the cascading plasma to order O {^^^ 
Equations of motion for the sound waves in the cascading plasma for generic ^ were 



derived in [29] . Previously, they have been discussed (solved) only to order O ( ^ 
[29]. Here, we extend the analysis to n = 4 in (2.13) at order 0(q°). 
Equations of motion for {-z^q ' -^/o ' -^i"o ' -^i^o ' -^k 0} form^° 



= 



3a;2 - 6a; + 2 
{x- 1)(,t2-2,t + 2) 



2x 



^J?,0 



+ 
+ 



32 



(n) 

8/3l,n 



16a;(2 - x) , 

K. 



16x^(2 -x)^ 



3(1 - x) 



2n 



{x-l){x^-2x + 2) 



^2 



2n 



(2.62) 



x'^-2x + 2 



+ 



X 



(n) , 3(X- 1) 



x^{2-xy 10a;(2-a;) 



+ 



x2(2 - x)2 ^^'O ^ 3x(2 - x)(l - x)2 + ' 



[2n] 



{nj 



+ 



+ 
1 



X 



Jn) 



(x-1) 



15(1 -X 



,3 ^'2n 



x2(2-x)2 '^'^ ' 5x{2-x) 

+ Jf^ , 



(n) 
^K,0 



in) 
^K,0 



(2.63) 



(2.64) 



10x(2 -x)(l -x)2 



i°We used (2.21) and (2.22). 
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in) 
%,0 



+ 



x-1 



(n) 
%,0 



+ 



+ 



X 



2{x - 1) 

x{2 - x) 



(n) 
^K,0 



(n) 



+ 



+ 

2 



3{x - 1)3 ^ 



3(x - 1) 



(2.65) 

(2.66) 



where the source terms J^^ , jI';^^} are functionals of the lower 

order solutions: z^^^, z''^\ z^J^fj , 4™o ' '^^m, 6m, V2m, Mm, C2m and /3i,„i, with m < 
n. Explicit expressions for the source term functionals are available from the author 
upon request. Apart from n = [34] and for {^^ q , z~Hfi 1 '^^.o) [29] these equations must 
be solved numerically. We use the same numerical approach as outlined in section 2.1. 



2.2.1 Order n^Q 
We find: 



,(0) 



2x — x^ ^ 



,(0) 



_ (0) _ 



(0) 
*,0 



,(0) 



0, 



(2.67) 
(2.68) 



2.2.2 Order n^l 
We find: 



.(1) 



^(1) _ 3:(2 - x) 
12(1 - x)^ 
2 
3 



ln(2a; — x'^) , 



Near the boundary, a; — > 0+, we have 

(1) _ _ 1 ,2/4,0 



— In a; 
60 



0{xHnx) , 



-^^x + x'/'sZ + 0{x'). 



Near the horizon, y ^ 0+, we have 



,(1) 



Altogether at this order we have 4 integration constants 

r 4,0 3,0 1,0 1,0 1 



which is precisely what is needed to specify a unique solution for {z^^l, z^j^} 



(2.69) 
(2.70) 

(2.71) 
(2.72) 

(2.73) 
(2.74) 

(2.75) 
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2.2.3 Order n ^2 

Near the boundary, x — > 0+, we have 

.(2) 



(2) 

Zfr, = X 



720 



(.^,-llna;)+0(x^/^), (2.76) 
Zh!o = ^ s^!^ + In x) , (2.77) 

1 1 o I ''^ 2,0 1 1 \ 2 3 , 2 / 4,0 , / 1^9 

m2 H mx ] x ' s, ^ + x [ s/n + 

480 40 ^'^ 32 J 15 '^'^ V V4320 

In 2 - -s'^\ - -4° - lCi'° - -4'?1 Inx + -^In^a;^ + 0{x'/') , 
20 ^'^ 60 ^ 72^^ 2 ■'^'7 720 J ^ ^' 

(2.78) 



= ^ f -T^sZ + i - i In 2 ) + x^/^ + 0{x' Inx) , (2.79) 



A 2,0 ^ 1 _ 1 
15^^'2+ 45 45 

4'o = ^(4',2+ (-^ + ^ln2 + 4^^1na;^ H-0(a;Mna;). (2.80) 
Near the horizon, y ^ 0+, we have 

4!o = 4\ + 0{y') , (2.81) 

4!o = (-2/3i,2 + 3204J - 84:; - I - + - f + l{X%r + Sei. 

(2.82) 

45 = 45 + ^(^') ' (2-83) 

41 = sZ + 0{y') , (2.84) 

4i = + ^iy") ■ (2-85) 

Altogether at this order we have 10 integration constants 

1-2,0 2,0 4,0 3,0 2,0 2,0 o 2,0 2,0 2,0-. (0 Ra\ 

\^K,2 1 ^H,2 1 */,2 ' *a),2 i *$,2 i °K,h ' , ^f^h ' °u,h ' ' \^.0\Jj 

which is precisely what is needed to specify a unique solution for {z^^oi 4qi 4o-> 
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2.2.4 Order n = 3 

Near the boundary, x — > 0+, we have 

41 = ^ (4> - ^ In 2 - 4;^) In + Oix'/^) , (2.87) 
4]o-xsZ + 0{x']n'x), (2.88) 



11 7 2 / 83 1 1 2,0 \ , 1 , ^ 

mi Sj^r, + m2 St^o ma; m 

180 80 V576 24 16 ^'7 64 



V9 15 ' 15 ' 15 "''^ J V V 360 

_ I ^4,0 2,0 I 239 2,0 301 2,0 ^ 2,0 , ^ !„ o _ !„ 4,0 

120^''' ^ 24 ^^'1 + 360 + 2160 2 + 4320^^ 120^-2 + 960 2^f'^ 

_ 1 2,0 _ 1 a2,0 _ 1 2,0 _ 1 2,0 _ 1 /•2,0 2,0 _ 1 2,0 i 9 _ 1 2,0 2,0 

2Q^K,3 72^2 24^*'2 60"' 12^^ ^^'2 45 10 ^^'^ 

In 2 st'? In 2 In^ 2 ^— ^ — | In a; + | In 2 H s'i^r. 

4 ■'^'^ 40 1440 103680 y V^SO 120 ^'^ 

+ J_^2,0 J_^2,o _ 1^2 \ 0(xV2i„^) , 

720^^ 360 2 17280y / ^ ^' 

(2.89) 

J3) 1 ^ 89 2 2,0 1 2,0 2 2,0 1 2,0 ,47 1 2p 

- rT5^^'^ ~ 405 + 9 '^'2 " 90'^'2 " 15 '^'3 " Yh'^^' + 540 ^ " 90 ^ 
1 
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lln2 4;^^+a;^/2 4'° + ^(x^ Inx) , 



(3) / 2,0 , / 1 o 5 2.0 , 1 2,0 I 2,0 1 2,0 I 1 i 2 o , „ 

4,0 = ^ ( «$,3 + ( 2/^1.2 - -Sk,2 + Y^Si:},2 + Sk,z + 2**,2 + 72 

+ J In 2 + In ^ + C'l^^' x) . 
2 ' 54 / / 



(2.90) 
(2.91) 



Near the horizon, y — > 0+, we have 

4\ = 4!^ + 0{y') , (2.92) 

4'^o = f |/9i,2 - 2/3i,3 + • • • V' + ^(^') ' (2-93) 
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where • • • denote dependence on lower order coefficients, except for {/9i,2, /9i,3} — tlie 
expression is too long to be presented here, 

45S = 45 + ^^(?/') ' (2-94) 



z, 



(3) _ ^3,0 



+ 0{y') , (2.95) 
= 4'°. + 0{y^) . (2.96) 
Altogether at this order we have 10 integration constants 

1-2,0 2,0 4,0 3,0 2,0 3,0 n 3,0 3,0 3,0 ; Qn\ 

which is precisely what is needed to specify a unique solution for {-z^Q' ^i?0' ^/^O' ^ufo' 

2.2.5 Order A 

Near the boundary, x — > 0+, we have 

1 on 265 



^(4) _^/^2,0 , /^_l/0 _ 1 In2 9-L'^e2,0 _ 1 2,0 _ 2,0 _-L 2,0 _ 

■^K,o — •^\^K,4'^ \ 2^'^'^ 12 ^^4'^i^,2 12 -^'2 2 144 

1 7 \ \ 

- -In2 4° + — ln2 1 lna; + — In^x ) + 0(x^/'^\nx) , 

2 24 y 48 y ^ ^' 

^g^o = x + C>(a;2 In^ x) . (2.99) 
Near the horizon, y — > 0+, we have 

4]o-4'H + 0{y'), (2.100) 
4!o = (^A,3 - 2A,4 + (-f - 484'° A^'° - f - 6404J - |a^'° 

+ mm'y - Ic^'' - ^4;° + 164:° + le^^. + 32^.^'^ + fixi'yy,,. (2.101) 

where ■ ■ ■ denote dependence on lower order coefficients, except for {/9i^2; /5i,3, /3i,4} — 
the expression is too long to be presented here. 

Altogether at this order we have 4 integration constants 

{4°4,4!!i,4^,/3i,4}, (2.102) 
which is precisely what is needed to specify a unique solution for {-z^q, -z^q}- 
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n 


1 
1 


9 


Q 
O 




2,0 
^K,n 




0.07891997 


-0.53177623 


1.48077259 


2,0 
^H,n 




-0.76150758 


2.66710077 


-5.15078326 


4,0 
*/,n 


-0.01641302 


0.02836594 


-0.00340767 






0.04361787 


-0.05326321 


-0.05682864 




2,0 




-0.03656488 


-0.55146784 




n,0 
^K,h 




0.14236466 


-0.61460152 


1.49213973 


n 

Pl,n 




0.33333333 


0.340767665 


-0.81625254 


n,0 


-0.00362335 


0.038494905 


-0.14888558 




n,0 


-0.00413161 


0.01711156 


-0.04003014 




n,0 




0.33416570 


-0.74238359 





Table 2: Coefficients of the normafizable modes of the sound quasinormal modes to 
order C»(qO). See (2.75), (2.86), (2.97) and (2.102). 



2.2.6 Integration constants for the sound quasinormal modes at C(q°) 

Here we tabulate (sec tabic 2) the integration constants for the normalizable modes of 
{^Kfi^ ^i?,o> ^/o ' ^i"o; -^i^o}- with n = {1, 2, 3, 4} obtained from solving the correspond- 
ing boundary value problems. 



2.3 Bulk viscosity of the cascading plasma to order O 

Equations of motion for the sound waves in the cascading plasma for generic ^ were 
derived in [29]. Previously, they have been discussed (solved) only to order O (^^^, 
[29]. Here, we extend the analysis to n = 4 in (2.13) at order 0{q^). 
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Equations of motion for {-sj^^ , z^^l , zj^l , z^^^ , z^^^^} take form^^ 



= 



(n) 



_Qx^2 
" (x^ -2x + 2)(x-l) 
2x^2 - xy 
3V2(a;-l)(a;2-2x + 2)2 



(n) 



+ 



X^ 



(n)' 
%,0 



+ 



2x + 2 
Ax{x - 2) 



32 



x^ 



2x + 2 



(n) 



31/2(3.2 _ + 2)2 -f^'O 



16(a;2-2a; + 4)(2-a;)V , lQ{x - 2)x , 

3V2(x2-2a; + 2)2^-^'0 ^ 3V2(a; - l)(x2 - 2x + 2)2 ^^n + 33/2(^ _ 1) '^2n 



, 64(2x2 -4a; + 3) („) 



B/?2,r 



8(3;2-2x + 4)/ji,„ 

31/2 (^2 _ 2a; + 2) 3^2 [x^ -2x + 2f 



rr\2n\ 



(2.103) 



(n) 
4,1 



+ 



X- 1 



(n) 
4,1 



8 (n) _ 3(a; - 1) 
a;2(2-a;)2 -^'^ 10a;(a; - 2) 



(n) 



+ 



(n) 
'/,0 



+ 



x2(2-x)2 -^'1 3V2(a;-l) 

^ ^ , q-\M 

3V2a;(l-a;)2(a;-2)'''"^^«'/ ' 



+ 



31/2 



10x(x - 2) 



(n) 



(2.104) 



= 



1 



X 



An) 



a; 



1) 



+ 



3V2(x - 1) 
31/2 



In) 



+ 



a;2(2-a;)2 "''^ ^x{x-2) 
3^/^ (n) 3V2 



(n) 



15x(a;-2)"^'0^ 15(x-l)3''2n 

A2n + J'ifc!^^ , 



10x(l - xY{x - 2) 



(2.105) 



(n) 



+ 



X- 1 



(n) 



+ 



2(x- 1) 
x{x - 2) 
2 



(n) 



31/2 (x - 1) 



(n) 



3V2:r(a; - 2) -^'O 3V2(x _ 1)3 ''2n + ^a,* > 



+ 



1 



X 



(n) 
^i^,l 



3V2(a; - 1) 



3V2(a; - 1) 



3"'2n 



, ^[2n] 



(2.106) 



(2.107) 



where the source terms { J"]^^', J'j^j', J^f"', ^7^^^', J'J^^-'} are functionals of the lower 

order solutions: 2;^-^, i -^$,1 ' ^k\i ^h})i "^7,0^' -^1,0^' '^l,d'' ^k})i ^Irai ^2rm V2rm 

^2m, C2m ^-nd {Pi,m, P2,m} , with ui < Ti. Explicit expressions for the source term 
functionals are available from the author upon request. Apart from n — [34] and 
for {z^^^i , Zff'-^} [29] these equations must be solved numerically. We use the same 
numerical approach as outlined in section 2.1. 



"We used (2.21), (2.22) and (2.67), (2.68). 
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2.3.1 Order n = 
We find: 



JO) _ (0) _ (0) _ (0) _ (0) _ ^ 

ho = 1 ■ 



2.3.2 Order n 
We find: 



^(1) _ ^(1) _ 



Near the boundary, x — > 0+, we have 



45 

= ^ - ^3^/2 Inx^ + Inx) . 
Near the horizon, y — > 0+, we have 

Altogether at this order we have 6 integration constants 

Sn"^'^ «2,o 1,0 1,0 1,0 1 (0 

which is precisely what is needed to specify a unique solution for {z^pi^ ^^i, ^^i}- 
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2.3.3 Order n ^2 



Near the boundary, x — > 0+, we have 



X (^0^2 + + C»(x3/2) , (2.119) 



,(2) 
-'HA 



X a 



l% + 0{x^\nx), (2.120) 



^(2) _ ™ J_q1/2 , ^ 2,0 J_ol/2 , 2 + J_„2,0 , J_ol/2 , \ _ ^ 3/2 3,0 
+ ^2 1^ 4,0 1 4 _ 1 2^0 _ I3I/2 _ 1^2,0 JL3I/2 In 2 + l3V\^° 

V V 2 •'^'^ 20 ^'^ 40 24 720 72 

(2.121) 

(2) ^ ^ f--a^'^2 - -3'/' - -4'° + -3'/' In 2^ + x'/' a'f, + 0{x' Inx) , (2.122) 



(2) _ / 2,0 , 



+ l^"^ + ^4,°i - ^3^^2 1112^ Inx) + ^(x^lna;) . (2.123) 
Near the horizon, y — > 0+, we have 

4\-atn + 0{y'), (2.124) 
- (-8a^. + f 3V24J _ 831/2^1,0 ^ ^3V2^o ^ _ ^31/2^0^ + I73V2 
+ 320 3^24.0 ^0 ^ + 24 31/2,1,0 ^0 ^ _ ^31/2^1 ^ _ 8 3i/2^2,c, ^ |^i,o _ 8,1.0 
- 2af^, + I3V2/J,, + 231/2^^^^ _ 1^3V2(^o^)2 _ ^3V2(A0 + 320a}J rj% 

(2.125) 

45 = 4;° + O(y^) , (2.126) 
4l = aZ + Oiy'), (2-127) 
4a = 4°. + (^(1/') • (2.128) 
Altogether at this order we have 10 integration constants 

/«2,o ^2,0 4,0 3,0 2,0 2,0 o 2,0 2,0 2,0-. roioQ^ 

which is precisely what is needed to specify a unique solution for -zj^^, zj\., zj^l, 
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2.3.4 Order n = 3 

Near the boundary, x — > 0+, we have 

(3) _ 
^K,l — ^ 



aZ + (-f 3^'^ + - aZ - ^^-) + 0{x^'') , (2.130) 

Zh\ = 2; a|°j + 0{x'^ In^ x) , (2.131) 

,(3) _ ^31/2^2^0 _ ^31/2 In 2 + 1 In 2 a^°2 

•^'^ V1O8O 160 ^'^ 160 ^'^ 360 80 ^'^ 

+ J_3i/2 2 + 14,0 + 1 In 2 aZ + -af, + -3^2/51 2 + -3^/2/52 2 - -af^ 
320 40 80 40 ^'^ 80 ^ ' 80 ^ ' 12 ^'^ 

_ J_4.o ^_Zl3i/2 + l3i/2in2 ^ lr4°2 - -4°il + -3V2ln2^^ 
240 V 576 32 16 ^^'^ 16 *'7 64 / 

1 3,0 1 3,0 1 r, 2 30 1 3,0 1 \ I 2/ 4,0 I / 3 1/2/3 

9«J,i - In 2 - — aj_2 - — aj_i Inxj + x + (^—3 / /?i,2 



360 ^'^ 360 ^ ' 720 1440 360 

_ il3 V2«:2,o J_3i/2^2,o 2 + J_3i/2 2C^° + 211^3^2 + -3 
432 ^ 120 ^ 144 518400 72 

60 ^ 1440 1440 360 ^'^ 120 ^'^ 40^^ 2 ■'^'^ 

^ 2,0 L/'^.O 2,0 ^ 2,0 2,0 ^2,0 ^ 1 r) 2,0 ^ 2,0 2,0 

24'^4>,2 -|^2^i '^^'2 12'^^ '^^'-'^ 20'^''^'^ 48 10'*^ '^^'^ 

(2.132) 

\ 90 15^'2 15 ^'^ 30 405 + 9^'^ ^ 90 

+ I3V240 _ 1 In 2 + -3^/2^1 2 + -3V2/?2 2 - —3^/2 1^ 2 + ^3^/2 1^2 2^ 
90 ^'^ 15 ^'^ 15 ^ ' 45 ^ ' 135 180 / 

+ x^/'' a^jl + 0{xHnx), 

(2.133) 



1 
360 
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(3) / 2,0 I / 1 q1/2 2,0 I 1 2,0 I 1 1 o 2,0 I 1 1 o 2,0 I 1 2,0 , 2,0 



°i - -3^/ 2 - -3^/V2 2 + -3'/' In 2 - i-3V2 In^ 2^ In x] 



_ Ilqi/2 _ ;^^2,o _ ^ 

54 3^'2 12^*' 

(2.134) 

Near the horizon, |/ — > 0+, we have 

4\-af,H + 0{y'), (2.135) 



{-ls'^'P2,2 + ^3V^/3,,3 + • • • + O(y^) , (2.136) 

where ■ ■ ■ denote dependence on lower order coefficients, except for {/32,2, /32,3} — the 
expression is too long to be presented here. 



zf}^af, + Oiy'), (2.137) 



zSl = al;l + 0{y'), (2.138) 
4i^aX, + 0{y'). (2.139) 
Altogether at this order we have 10 integration constants 

r 2,0 2,0 4,0 3,0 2,0 3,0 p 3,0 3,0 3,0-. M -\ 4n\ 

\«i^,3 ' (^H,3 ' «/,3 ' «a;,3 > %,3 ' (^K,h ' ^2,3 , df ^ , d^^h ' «$,?iJ ' (^^.i4Uj 

which is precisely what is needed to specify a unique solution for {z^^\, z^^\, ^j?], z^\, 

2.3.5 Ordern = 4 

Near the boundary, x — > 0+, we have 

4\ - - {aZ + (-^4!^ - - aZ - §3^/^ In 2 + ^a^^^ + + 1^^%,, 

+ ^3 In^ 2 + g3 + L:,y^sZ - ^ In 2 a^, - 1 In 2 a^'° + 13^^/3,,) In x 

- +0{x^/Hnx), 
48 J 

(2.141) 

= a; a^°j + 0{x^ In^ x) . (2. 142) 
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Near the horizon, y ^ 0+, we have 

4!i = 4> (^(?/') , (2.143) 

+ 16 3V^.i° + H^3V^4J + f 3VV,. + |3V^A?,, + f 3V^)/3., + • • 
+ 0(1/') , 

(2.144) 

where ■ ■ ■ denote dependence on lower order coefficients, except for {^2,2, (^2,3, ^2,4} — 
the expression is too long to be presented here. 

Altogether at this order we have 4 integration constants 

r 2,0 2,0 4,0/3-1 1 A 

{«i,4 > «//,4 > > /^2,4} , (2.145) 

which is precisely what is needed to specify a unique solution for {z^^\, z^^\}. 

2.3.6 Integration constants for the sound quasinormal modes at 0{(\^) 

Here we tabulate (see table 3) the integration constants for the normahzable modes of 
{"^^Kii ^1/1' -^/"i ' -^i"!}- "^^^ ^ = {1) 2, 3, 4} obtained from solving the correspond- 
ing boundary value problems. 



3 Challenges of computing transport coefficients to all orders 




In the previous section we detailed the computation of the speed of sound and the bulk 
viscosity of the cascading plasma, perturbatively in ^ (note that ^ ~ (in ^) ^ for^^ 
T ^ A (1.4). ) The results of the analysis presented in section 4 indicate that although 
^ (^) psrturbative expansion is in excellent agreement with the full (nonperturbative 
in ^ computation for in the high temperature regime, this expansion does not 
converge below T ~ (1---1.5)A, which is about twice as high as the temperature of 

^^For exact temperature dependence of K^, see [22]. 
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1 
1 


9 
z 


Q 
O 




2,0 




0.01069638 


0.23174494 


-0.98113866 


2,0 
«i?,n 




0.08509403 


-0.78349123 


3.54269325 


4,0 


0.03737595 


-0.06797689 


0.08231942 




n 


-0.08394261 


0.10167870 


-0.05749855 




2,0 
%,n 


-0.29631916 


0.38704391 


-0.14363558 




n,0 




-0.18958164 


0.69624351 


-1.81005605 


n 

P2,n 




0.13225837 


-1.69770959 


2.26988336 


n.O 
^f,h 


0.00580530 


-0.05087126 


0.17559889 




n.O 


0.00630707 


-0.02271965 


0.05689077 




n,0 


0.11330816 


-0.38874716 


0.95684007 





Table 3: Coefficients of tfie normalizable modes of the sound quasinormal modes to 
order C»(qi). See (2.118), (2.129), (2.140) and (2.145). 

the deconfinement phase transition. Thus, using perturbative analysis only we can 
not compute the bulk viscosity of the cascading plasma at the transition point. In 
this section we would like to explain the difficulty in going beyond the perturbative 
analysis as the latter might effect the analysis of other quasinormal modes in the 
cascading plasma, specifically those that could be responsible for the chiral symmetry 
breaking transition [24] . 

To understand the problem, it is instructive to go back to the numerical com- 
putation of the cascading plasma equilibrium equation of state. This was solved both 
perturbatively and non-perturbatively in ^ in [22] . On the dual gravitational side this 
computation involves finding the black hole solution in asymptotic KT geometry, i.e., 
determining the gravitational fields {h , f2 , fs , K , g} (2.3)-(2.7). Above gravitational 
fields have non-normalizable (in some cases singular) modes both near the horizon 
X — > 1_ and near the boundary x ^ 0+. Thus numerical integration must be done on 
an open interval x e (0, 1), i.e., we need to provide the boundary conditions for the 
gravitational fields as the series expansion in x near the boundary and m y = 1 — x 
near the horizon. These series expansions must be fairly precise since, for examples 
the coefficient of the normalizable mode for /2 (dual to the vev of the dimension eight 
gauge invariant operator of the cascading plasma) enters at order x^ near the bound- 
ary, which is subdominant to coefficients x'*/^ In'^(x) with n = 0, 1,2,3 of the general 
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boundary expansion [22] 

oo n 

f2^ao + J2Yl ^ ■ (3.1) 

n=l fc=l 

What saves the day, and ultimately allows for the full non-perturbative computations, 
is the fact that at each fixed order n, the maximum power of Inx in (3.1) happens to 
be bounded, k < n. Thus, the series expansions of the type (3.1) are just generalized 
Taylor series expansions, which can be easily determined to any given order in n — 
the total number of expansion coefficients at order n grows as 0{n'^). The situation 
would have been completely different, had the summation of k extend to infinity. Here 
one would have to solve exactly for the series in (P^ In x) at each order in x. Given the 
complexity of the equations involved the latter appears to be impossible. 

Unfortunately, precisely this problem occurs in computation of the quasinormal 
modes in the sound channel. Consider for example the gauge invariant fluctuation Zf 
(2.8). As for the gravitational field /2, it depends on the vev of the dimension eight 
operator — so its exact boundary asymptotic can not be specified with an accuracy of 
less than 0{x'^Y^. Collecting (2.71), (2.78), (2.89) we find 

^^'°=<-80^^"^+ 720^^"^-64K!^^^ + -- j+^i^^^^ 

(3.2) 

where we explicitly indicated only the leading In a; dependence at each order We 
further verified that z^^^ — O (^x In^ x) as a; — 0+, and that in fact all the perturbative 
expansions for {zh, Zf, z^^, z^, zk} do not truncate in k — for example, 

oo oo 
n=2 k=0 

It would be interesting to develop computational techniques to deal with this difficulty. 
Notice that the high temperature perturbative expansion provides an effective cutoff 
on the power of In x in the boundary asymptotics since each additional factor of In x 
comes with a factor of 



^^Of course, in order to get reliable numerical results boundary asymptotics must be more precise 
— in our high temperature analysis we used expansions to order 0(a;^/^), which is five more orders 
beyond the highest order at which the normalizable coefiicients of the fiuctuations enter. 
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4 Perturbative hydrodynamics of the cascading plasma 



In this section we present results of the perturbative high temperature analysis of the 
speed of sound waves and the bulk viscosity in the cascading plasma. We begin with 
discussion of the consistency checks on our analysis. Next, we move towards discussion 
of comparison between exact speed of sound (as given by (2.60)) and its perturbative 
high temperature expansion. This will allow us to comment on the convergence prop- 
erties of the high temperature expansion. Extending the numerical analysis of [22] we 
show that chirally symmetric deconfined phase of the cascading gauge theory plasma 
becomes perturbatively unstable below the critical temperature of the deconfinement 
transition TanstaUe = 0-87A9{0)Tcriticai- We comment on the possible source of the in- 
stability. Finally, we discuss the bulk viscosity bound of [31] for the cascading plasma. 

4.1 Consistency of analysis 

4.1.1 The first law of thermodynamics 

Cascading gauge theory plasma has a single scale A. It only makes sense to discuss 
the thermodynamics/hydrodynamics of the theory, provided one keeps A fixed. As 
explained in [22], enforcing that A is temperature independent leads to the following 
condition 

^(i^.-^lna^)=0, ao^^lim/.. (4.1) 

From (2.4) we find 

1 P2 / 1 _ 1 \ / 1 , o „ 49 17 



an =do ( 1 - h( — ln2-— — In2 2-F^^^ ^ln2, , 

° °' V18 12) V54 648 216 J 



■k 



+ 7 + 77777 In 2+ ln^2 \n^2] — + 



64 243 162 216 J \ 



■k 



(4.2) 



Next, we compute the temperature of the black hole (2.1) using (2.3)-(2.4), and an 
explicit expression for Gxx — given by eq. (2.6) of [22]. We further invert the temper- 
ature relation to obtain oq, and ultimately from (4.2) the perturbative expression for 
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+ (I - |ln2 - 3e°, - f «,)^ - ^(A?,,)^ - 6^3% + ^ln^2 - ^A?,, A°, 



(4.3) 



Using (4.3) we find from (4.1) 

T 

dA=0 



dlnT 



pA p6 p8 /P^0~ 



+ 



p 



10 



(4.4) 



Given (2.57), (2.58) and the expression for the entropy density s of the black hole 
(2.1), the first law of thermodynamics 



dV = sdT , 
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n 


1 


2 


3 


4 









8.2 X 10-*^ 


-3.5 X 10-^ 


-6.1 X 10"^ 



Table 4: Comparison between and of (4.4) and (4.5). 



leads to an alternative expression for 



d\nT- 



dlnT 



pA p6 p8 / P^^ 

dV-sdT=0 



-2P' + (-^ + 2c^° + 4.f + in2) g + (4.f - l^r - - Ic.^' + g 

-60(<J^.2C- + 24^1n2.iln^2-31n2.1n2C--24eo,)|! 

- 7«f In 2 - ^ In^ 2 + In 2 C' '° + ^ 2 Ci '° + 2 In 2 «f + In^ 2 «f + ^ In' 2 
+ 40(r,°J^ + ^(A?,J^ - 36^3% - 180<, v% + f (A?,J^ - f A?,, A^,, 

- 240«,)=' + |(A° J' + 180^2% - 301n2 (^J^ - 121n2 ^^^^ - ^ ln2 (A°,)^ 

- ^ - 120«J^ - 60«J^ - 48e,% - 24e,% - f (A?,J^ 

p8 /P^*'' 



(4.5) 



Comparison between ^^d provides a highly nontrivial test on the consistency 
of the analysis. Using the data from the table 1, the results of such comparison are 
presented in table 4. 



4.1.2 (?g from the equilibrium thermodynamics and from the hydrodynamics 

Our second consistency test compares the predictions of the equilibrium thermody- 
namics for f3i^n — Pf^n^"^" from (2.61) with the direct computation of = /'ij*""""' (see 
(2.14) for the parametrization and table 2 for the results) . Using the data from the 
table 2, the results of such comparison are presented in table 5. 
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n 


1 


2 


3 


4 




( lothermo \ 

1 ' 1 
\ ^1." / 





-1.8 X 10-^ 


4.0 X 10-^ 


5.5 X 10-^ 



Table 5: Comparison of coefficients Pi^n in the sound wave dispersion relation (2.14) 
with the predicted values from the equilibrium thermodynamics, (2.61). 



4.2 Speed of sound and perturbative instability of deconfined chirally sym- 
metric phase of the cascading plasma at low temperatures 

The speed of the sound waves can be computed from the dispersion relation of the 
quasinormal modes in the sound channel. In the high temperature expansion it is 
given by (see (2.14)) 



If / 4\ p4 / 4^ \ P6 



high-T 



3K^ \ " ' 9/ V 3" 'V i^l 

+ (2fl,.-ifl,3 + (A.)^)g + oO} 



(4.6) 



where the coefficients are given in table 2. Alternatively, it can be evaluated from 
the equilibrium thermodynamics for any temperature (2.60), 



thermo 



dV _1 7- 12a,,o - 6P2 ^ 
dS 3 7 + 4a2,o + 2P2 ^ ■ 



Prom data in table 5, we see that there is an excellent agreement between (4.6) and the 
high temperature expansion of (2.60). Figure 1 represents comparison between (2.60) 
and (4.6) over a wide range of temperatures, i.e., not necessarily when ^ <^ 1. The 
blue dots represent the speed of sound computed from (2.60), slightly improving the 
analysis in [22]. The solid lines represent successive high temperature approximations 
to the speed of sound wave in the cascading plasma (4.6) to orders O (black), 

C (^^j (purple), O j (green) and O j (blue). It is convenient to plot the data 
with respect to = ^ — | In 2, rather than with respect to ^. The vertical red hne 
represents the deconfinement temperature of the cascading gauge theory plasma. No- 
tice that the second C and the higher orders of the high temperature expansions 
become indistinguishable with the exact numerical data (blue dots) for > 3. Using 
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Figure 1: (color online) Speed of sound (blue dots) vs. its successive high temperature 
approximations. The vertical red line represents the temperature of the deconfinement 
phase transition in the cascading plasma. 



results of [22] 



T T 

k, = 3 <^ v~l-43 fti2.34, 

A 7^ 

^ -'■ critical 

T T 
k, = 2 -^1.00 ^1.63, 

^ -'■ critical 



which suggests that the high temperature expansion converges for temperatures above 
T> (1---1.5)A. 

Numerical analysis of the equilibrium thermodynamics in [22] were done to tem- 
peratures only slightly below the deconfinement temperature. Here, we extend the 
computations to lower temperatures. Notice from figure 1 that the speed of sound 
squared Cg appears to cross zero (and turns negative) for fc^ < — 1. A more detailed 
analysis presented in figure 2 show that this is indeed so. A speed of sound vanishes at 

cl{kunstable) = =^ kunstahle = —1.230(3) . (4.9) 

While 

~ (^s ~ ^unstable) the viciuity of the instability, we find 



1/2 

Me \ 



2 I 1 ^unstable \ / a -\ r\\ 

c, ~ 1 — , (4.10) 
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Figure 2: Speed of sound in the vicinity of the perturbative instabihty of the cascading 
plasma. 



where 

^""^'""'^ = 0.53728(8) , or T^nstaUe = 0.8749(0)T,„w . (4.11) 

In fact, the near- unstable thermodynamics of the cascading theory is rather interesting. 
In figure 3 we show the free energy density as a function of ks (left plot) and as a 
function of ^ (right plot). Notice that there are two phases separated by a continuous 
phase transition at ks = kunstabie- The phase with kg > kunstabie (which is continuously 
connected to a high temperature deconfined chirally symmetric phase of the cascading 
plasma) has a lower free energy compare to a phase with ks < kunstabie — the latter two 
phases are degenerate in temperature (see the right plot on figure 2) with the limiting 
temperature being reached precisely at ks = kunstabie- It is clear from the right plot of 
figure 3 that^^ 

is finite. Thus, the only way the speed of sound can vanish at T = TunstaUe is if the rate 
of change of the energy density diverges at the unstable point. The plots in figure 4 
show that this is indeed the case. 

Vanishing of the speed of sound as in (4.10) implies that the specific heat cy of the 
cascading plasma diverges near the unstable point with the critical exponent a = 0.5^^: 



s 

Cv = ^ (X 

ci 



Tunstable 



T 



-1/2 

(4.12) 



^"^We did further, more detailed, numerical analysis to confirm this. 

-"^^This coincides with the mean- field critical exponent a at the tricritical point [35]. 
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Figure 3: (color online) The free energy density in the vicinity of the perturbative 
instability of the cascading plasma. The red dots on the right plot correspond to 
< Knstabie and the blue dots correspond to kg > KnstaUe- 

Exactly the same critical behavior was found in A/" = 2* plasma with mass deformation 
parameters nif < nif, [19,31]. 

Whenever < 0, the thermodynamic system is unstable with respect to density 
fluctuations. It is interesting to understand the source of this instability. Recall that 
the zero temperature the vacuum of the cascading gauge theory spontaneously breaks 
chiral symmetry. Thus, it is conceivable that the perturbative instability observed 
in the equilibrium thermodynamics of the deconfined chirally symmetric phase of the 
cascading plasma is associated with the formation of chiral condensates. We comment 
more on this in the conclusion and for further analysis refer to future work [24]. 

4.3 Bulk viscosity bound in the cascading plasma 

The primary goal of hydrodynamic analysis of the cascading plasma presented here 
was to verify the bulk viscosity bound in strongly coupled gauge theories conjectured 



where p is the dimension of space of the plasma. It is known that (4.13) is satisfied 
in all explicit realizations of gauge theory/string theory correspondence [19,31]. The 
bound is saturated for all Dp branes [31,36]. It is also known that one can engineer 
phenomenological models motivated by gauge/string correspondence that would violate 
the bulk viscosity bound [37]; finally, the bound is violated in weakly coupled gauge 



in [31]: 




(4.13) 
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Figure 4: (color online) The energy density in the vicinity of the perturbative instability 
of the cascading plasma. The red dots on the right plot correspond to kg < kunstaUe 
and the blue dots correspond to ks > kunstaUe- 



theory plasmas 

What makes cascading gauge theory plasma interesting in this context is that it 
provides an example in the framework of the gauge/string duality where the bound 
(4.13) is saturated to leading order in the high temperature expansion [29]. As we 
explained in section 3 our current computational framework is inadequate to compute 
bulk viscosity of the cascading plasma at low temperatures. In the high temperature 
expansion it is given by (see (2.14)) 

Using (4.6) is can be rewritten as 
C , , /^27^ , 3 , 27^ \ 2 /243^ 729^ ^ «t ^ 



729^ ^ 81^ 243^ V, A215 ^ 32805^ ^ 98415^ 3 



128 ' ■ 64 ' 128 '7 V1024 ' 2048 ' ' 4096 

6561 ^ 2 2187^ 135 2187^ 6561 ^ ^ 98415^ ^ 2 

H A 2 H A 3 H ^ /5i 4 H /52 2/^1 3 H /32 2/^1 2 

1024'^ ' 1024'^ 512 1024^ ' 1024^ ' '^ '^ 4096 ^^'^-^^'^ 

6561 1215 2187 19683 2187 \ 4 

+ Y5^/^2,2/^l,2 + + Y5^/:^2,4 + ^5^/^2,3/^1,2 + Y5^/^2,3 J 

+ O (5^) , 

(4.15) 
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0.1 0.2 0.3 0.4 0.5 0.6 



Figure 5: (color online) Successive high temperature approximations for the bulk vis- 
cosity bound for the cascading plasma. The viscosity bound implies ^ > f^- The solid 
vertical red line represents the temperature of the deconfinement phase transition in 
the cascading plasma. Vertical dashed lines indicated the expected convergence of the 
high temperature expansion: the left line corresponds to T ^ 2.34Tcriticai and the right 
line corresponds to T ^ 1.63TcriUcai- 



where we introduced 

5^2Q-c^). (4.16) 

The results of the high temperature computations are presented in figure 5. The solid 
lines represent successive high temperature approximations to the bulk viscosity bound 
in the cascading plasma (4.15) to orders O j (purple), O (^^^ (green) and O j 

(blue). Recall that the bound is exactly saturated at order O (^^^- The vertical solid 
red line represents the value of 5 = Scriucai at the deconfinement phase transition. The 
dashed red lines indicate the expected convergence of the high temperature expansion 
of the hydrodynamic quantities deduced from the high temperature expansion of the 
speed of sound, see figure 1. 

From figure. 5 we see that at least in the high temperature expansion the bulk 
viscosity bound (4.13) is satisfied. Since the deconfinement phase transition in the 
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cascading plasma is of the first order [22], we do not expect any singular behavior in 
the bulk viscosity [32] in the vicinity of the transition. Notice that there is almost no 
difference between O (y^s^ (green) and O (yj^ (blue) approximations to the viscos- 
ity bound all the way to the deconfinement temperature Tcruicai- This suggests that 
the high temperature expansion for the bulk viscosity might have better convergence 
properties than that of the speed of sound. If we take the high temperature results at 
the deconfinement transition seriously, we find that 



-0.6(1). (4.17) 

deconfinement 



QCD slightly above the deconfinement phase transition is nearly conformal. For cj. 
in the range 0.27 — 0.31, as in QCD at T = 1.5Tdeconfinement [39,40], we are well inside 
the expected validity range of the high temperature expansion, resulting in 



0.05 -0.14 (4.18) 

QGP 



for the cascading gauge theory plasma. 



5 Conclusion 

In this paper we presented detailed analysis of the transport properties of the deconfined 
chirally symmetric phase of the cascading plasma at strong coupling, using the gauge 
theory/string theory correspondence. We developed the high temperature expansion 
to order ~ i^^li) ^ ^'-'^ thermodynamic and the hydrodynamic properties of 
the theory and identified challenges in going beyond perturbative in Q hydrodynamics. 

We computed the high temperature expansion of the bulk viscosity of the cascading 
plasma. We showed that the bulk viscosity bound proposed in [31] is satisfied in such 
plasma. We argued that results for the bulk viscosity are hkely to be reliable up to the 
deconfinement temperature with bulk viscosity being about 60% of the shear viscosity 
right at the deconfinement transition. Much like in other holographic models of gauge 
theory/string theory duality [19] we observe a rapid drop in the bulk viscosity above 
the deconfinement transition. 

An interesting byproduct of our hydrodynamic analysis was the discovery of the 
perturbative instability of the deconfined chirally symmetric phase of the cascading 



39 



plasma. Specifically, extending analysis of [22] we identified a continuous phase transi- 
tion in the plasma (slightly below the critical temperature of the first order deconfine- 
ment transition) where the speed of sound squared vanishes and becomes negative. A 
similar phase transition was observed previously in A/" = 2* gauge theory plasma [19,31]. 
Although in the former case it is difficult to speculate as to the origin of the instability, 
it is tempting to relate the same instability in the cascading plasma with the develop- 
ment of the chiral condensates responsible for the breaking of chiral symmetry. The 
fluctuations of such condensates are massive at high temperatures [24]. Exactly for 
this reason there is no high temperature regime for the deconfined cascading plasma 
with broken chiral symmetry — correspondingly, there can not exist a black hole solu- 
tion on the warped deformed conifold with fiuxes [2] at high temperatures. Of course, 
this docs not exclude the possibility of such a black hole solution at low temperatures. 
The hydrodynamic stability of the symmetric phase all the way down to the decon- 
finement transition suggests though that the existence of the black hole in the broken 
phase would not modify the cosmological scenario proposed in [41]. We return to these 
questions in more details in future work [24] . 
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